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We investigate the performance of inhomogeneously broadened spin ensembles as quantum mem- 
ories under continuous dynamical decoupling. The role of the continuous driving field is two-fold: 
first, it decouples individual spins from magnetic noise; second and more important, it suppresses 
and reshapes the spectral inhomogeneity of spin ensembles. We show that a continuous driving field, 
which itself may also be inhomogeneous over the ensemble, can enhance the decay of the tails of 
the inhomogeneous broadening distribution considerably. This fact enables a spin ensemble based 
quantum memory to exploit the effect of cavity protection and achieve a much longer storage time. 
In particular, for a spin ensemble with a Lorentzian spectral distribution, our calculations demon- 
strate that continuous dynamical decoupling has the potential to improve its storage time by orders 
of magnitude for the state-of-art experimental parameters. 



I. INTRODUCTION 

Ensembles of a large number of emitters, e.g. cold 
atoms, polar molecules, and electronic spins, are very 
appealing systems for quantum memories [1] and quan- 
tum repeaters [2] for long distance quantum communi- 
cation. Recently, experimental progress has led to the 
demonstration that hybrid quantum systems composed 
of superconducting qubits and a collection of electronic 
spins represents a promising new route [3-7] , which takes 
the advantage of long coherence times in spin ensem- 
bles and fast operation of a superconducting qubit. The 
large number, N, of spins provide a collective enhance- 
ment of the interaction strength between the spin en- 
semble and the electromagnetic field, and the collective 
coupling strength scales as = (J]]^ QkY^'^ ^ VNq. The 
strong coupling regime for an ensemble of electronic spins 
and a superconducting resonator has been experimentally 
demonstrated recently [3-6], i.e. ^/Ng ^ 7, where 7 
and n are the spin's and resonator's damping rates. 

Despite these achievements, the realisation of a quan- 
tum memory for microwave photons using spin ensem- 
bles remains a challenge due to both, the linewidth of 
the resonator and of the spin ensemble [4]. Apart from 
the dephasing of individual spins due to magnetic noise, 
the spectral inhomogeneity over the ensemble will lead to 
the decay of the collective superradiant spin wave mode 
of the spin ensemble to dark states. The performance of 
a quantum memory is largely dependent on the proper- 
ties of the spectral inhomogeneity [8, 9]. In particular, if 
the tails of the spectral distribution decay faster than a 
Lorentzian profile, the cavity can provide coherence pro- 
tection and improve the storage time of a quantum mem- 
ory due to the collective enhancement of the coupling 
strength. However, Lorentzian distributions are found in 
many solid-state systems, e.g. electronic spin ensembles 
in diamond and impurities in crystals. In these systems, 
the memory time will still be limited significantly due 
to the spectral inhomogeneity and would hardly benefit 
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Figure 1. (Color online) Setup of a cylinder- like diamond 
sample coupled with a high-Q quantum resonator. The dy- 
namical decoupler drive resonator generates continuous mi- 
crowave driving field. The driving field at the position with 
a distance r to the center of the sample is calculated using 
the Biot-Savart law. The radius and height of the sample are 
denoted as d and h respectively. The distance from the center 
of the sample to the drive resonator is R. 



from the strong collective coupling [8, 9]. 

In this work, we investigate how continuous dynamical 
decoupling can be used to construct a memory based on 
driven spin ensembles that has a long storage time. It is 
known that continuous driving fields can decouple spins 
from magnetic noise [10-14]. How continuous dynami- 
cal decoupling can improve the performance of quantum 
memories has not yet been examined in detail. Here, 
we show that besides the protection against magnetic 
field noise, continuous driving can have another promi- 
nent role in the present context, namely it can suppress 
and reshape the spectral inhomogeneity of spin ensem- 
bles. For Lorentzian distributions, continuous driving 
leads to a much faster decay of the tails of the effec- 
tive spectral distributions. This effect makes all the spin 
transition frequencies far from the energy gap induced by 
the strong collective coupling, and thus the driven spin- 
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ensemble memory can benefit from the cavity protection 
[8, 9]. To analyze the performance of the driving sys- 
tem as a quantum memory, we calculate the decay of the 
collective Rabi oscillation and the state fidelity after a 
certain storage time by solving the Schrodinger equation 
with the Laplace transformation. Our numerical calcula- 
tions demonstrate that the storage time of a driven spin 
ensemble can be improved by orders of magnitude with 
reasonable experimental parameters, even if the continu- 
ous driving field itself is inhomogeneous over a finite size 
of sample. 



the following, our calculations will use the conservative 
estimated driving parameters: Bmin = ^^^^^"^ = lOA 



and 



2R 



= 5mm(l + 0.05). We 
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remark that the design of the '8' shape high Q quantum 
resonator, see Fig.l, is to reduce the effect of the driving 
fields on the cavity mode since the the total driving field 
on the resonator should vanish. 



III. EFFECTIVE SPECTRAL DISTRIBUTION 
OF DRIVEN SPIN ENSEMBLES 



II. SETUP OF DRIVEN SPIN ENSEMBLES 

We consider a cylinder-shaped diamond sample, which 
is coupled to a superconducting resonator, see Fig.l. The 
radius and height of the sample are denoted as d and h 
respectively. The amplitude of the microwave driving 
field at the point with the distance r to the center of the 
diamond sample, can be estimated using the Biot-Savart 
law B{r) = 9^mL J and we get 



B{r) 



Jo 



R{R -r COS 0) 



_2i^rcos6> + r2)3/2 



which can be written as 



Bir) = 



2R 



lie, (1) 
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where dl is the differential element vector of the drive 
resonator in the direction of current, fio is the magnetic 
susceptibility, R is the radius of the dynamical decou- 
pler drive resonator, and / is the amplitude of the cur- 
rent. It can be seen from Eq.(2) that when the condi- 
tion R > Ad is satisfied, the inhomogeneity of the driv- 
ing field amplitude will be less than 5% (i.e. the rela- 
tive difference between the maximum and the minimum 
values across the sample). For diamond samples, it is 
feasible to achieve a concentration of Nitrogen -Vacancy 
(NV) centers of po = lO^^mm"^ [15]. For a thin dia- 
mond film with the radius d = 0.1mm and the depth 
h = 0.1mm, the number of NV centers is around 10^^. 
It has been estimated from experiments that the inter- 
action between an individual spin and a superconduct- 
ing resonator is of the order of lOHz [4] and thus a col- 
lective coupling strength on the order of lOMHz can be 
achieved. Assuming A the value of the width of the spec- 
tral distribution (which is of the order of MHz [4]), we 
can estimate that the required driving current is about 
0.3 (A/MHz) A in order to generate a microwave field 
(which creates a Rabi frequency) of lOA on the diamond 
sample with a driving source of radius R = 0.5mm. The 
required current can be reduced by using several res- 
onator loops. These parameters are within reach of cur- 
rent experimental technology. Careful optimization may 
easily improve homogeneity and strengthen the drive. In 



As we have discussed in the above section, the mini- 
mum and maximum values of the amplitude of the contin- 
uous driving field on the diamond sample are B^iiji and 
Bmax respectively. We denote the probability function 
of the driving amplitude as f{B), and have f{B)dB = 
po • {27Trdr)h where po is the concentration of NV center 
and h is the depth of the diamond sample. We remark 
that the actual shape of the diamond sample will not 
affect the main conclusions of the following discussions. 
After some calculations, we can see that the distribution 
of B actually has a rectangle profile, namely 



f{B)dB = po{27Trh)dr = 



AnipoR^h) 

'^Bryj.jY)^ 



dB. 



(3) 
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Thus, the normalized probability function for B 
f{B) = 1/Ab for B e [Bmin, Bmax], whcrc Ab = 
Bmax — Bmin is the inhomogcueity of driving field. The 
local Hamiltonian of an individual electronic spin is 



^ (k) 
2 ' 



• Bkai^^ cos (ujot) 



(4) 



where ujo is chosen to coincide with the central value of 
the energy spectrum (which corresponds to the zero field 
splitting of NV center ground state namely 2.87GHz and 
shifted by the applied magnetic field) and ai^^ , ai^^ are 




Figure 2. Lorentzian spectral distribution (blue) and the one 
modified by continuous driving (red). We use the parameters 
Bmin = lOA and Bmax = 10. 5A. In the red curve the zero of 
frequency is located at {Bmin + Bmax)/ 2. 
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Pauli operators. In the interaction picture, the effective 
Hamiltonian is written as 



Afc 



.(fe) 



(5) 



Taking into account the distribution of A/., one can show 
that the correction of 4^ is neghgible when ^ 
[16], and thus the effective couphng strength can be ap- 
proximated as gk = (i.e. Q = C^j^giy^'^ = ^/2). 



where we denote = ojk — ojq. After diagonahzing the 
above effective Hamiltonian, we obtain a dressed two- 
level system with 



IV. TRANSMISSION FUNCTION AND 
COLLECTIVE RABI OSCILLATION 



\i)k{i\) 



(6) 



with ujk = \/Al^Bl > Bmin where W)^ and \\)^ de- 
note the effective spin up and down states (namely the 
eigenstates of the effective Hamiltonian in Eq.(5)). The 
value of Bj. is uniformly distributed within the range 
[Bmin^Bmax\i and wc cousidcr a Lorentzian distribution 
for A/e, namely 



p(Afe) 



2A 



^(A2+4A|)- 



(7) 



Thus, we can calculate the normalized probability distri- 
bution of the effective energy spectrum as follows 



= 2 j KB) 



Ak = 



S2 dAk 



dB 

(8) 



which leads to 



A2 + 4(^2_;^^Srf^fc 



After performing the integration, we get 



(9) 
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^ Ok e {Bmax, OO) 

(10) 

where ^{B,Ok) = t^ii-\B A / ^ {Auol + A^XcD^ - B^)) 
and Ab is the inhomogeneity of the driving field. 
In the ideal case of perfectly homogeneous driving, 
namely when Bmax = Bmin^ th e distribution function is 
p(cDfc) = 4u,,A/7T^Lol - Bl^^ [A2 + 4 _ ^2^^ 
UJk > Bmin- The interaction between individual spins 
and the field for a resonant cavity is 



(ai^^ ± iG^^)/2 and gk ^ 



(11) 



where a^l^ = ((7^' ± ^cr^"0/2 and gk - lOHz. With 
the rotating wave approximation, which is valid when 
gk <C Xo, we can obtain the effective interaction Hk-c = 
^/c(<^ lt)fe Ul li)/c (tl) between the cavity field and the 
dressed two-level system with 
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gk 



z9k 1 



Ak 

Ok 
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We consider an ensemble of electronic spins coupled 
to a resonator, whose frequency and damping rate are 
denoted uOc and n respectively. When the number of ex- 
citations is low as compared to the total number of elec- 
tronic spins, the spins can be described as a collection of 
oscillators in the Holstein-Primakoff approximation with 

^ and cr^^ = 



total Hamiltonian can be written as 



a\ak 



at, thus the 



AT 



N 



H = ujca\ac + ^cj/ctt^a/e + ^^gk{a\ac + aka\) (13) 



k=l 



k=l 



Under the driving of continuous field, we have a collection 
of dressed two- level systems coupled with the resonator, 
and the effective Hamiltonian is 

AT AT 

H = uJcdlac -^^uJkalak + ^ y (4^c + ^/c4) (14) 

/e=l k=l 

The difference between Eq.(13) and Eq.(14) is that the 
original Lorentzian distribution of uJk is now modified to 
the new distribution of cD/e, as shown Fig. 2. The width 




Figure 3. Transmission of a resonator (as in Eq.(15)) coupled 
to a Lorentzian spectral distribution of spins without (a) and 
with (b) continuous driving. The parameters are Bmin = 
lOA and Bmax = 10. 5A, and k = O.lMHz, 7 = 10"^MHz. 
(c)The collective Rabi oscillation between the cavity and the 
superradiant spin wave mode with the initial state |1, G). The 
curves show the quantity F(t) = | (1, exp (—itH) |1, C) |^ as 
a function of time. 
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of the distribution depends on the inhomogeneity of the 
driving field, and most important the decay of its tail 
becomes much faster. This fact will greatly improve the 
performance of such a quantum memory as we will show 
in the rest of this work, although the effective coupling 
is decreased by half. 

We assume that the cavity is on resonance with the 
spins, namely its frequency is uOc = for the case with no 
continuous driving; while the cavity resonant frequency is 
cjc = ^0 + {Bmin-^Bmax)/'^ Corresponding to the dressed 
spin wave mode under continuous driving. The transmis- 
sion function of the cavity from the standard input-out 
theory is given by [9, 17] 



IK 



2 (ujQ — iti/2 — uj) 



k 



Z7/2 ■ 



(15) 

where 7 is the damping rate of spins. In the strong 
coupling regime, namely the collective coupling is much 
larger than the decay rates (i.e. Vt ^ 7,/^), the exci- 
tation is coherently exchanged between the cavity and 
the superradiant mode of the ensemble of emitters. The 
so-called vacuum Rabi splitting corresponds to the ap- 
pearance of double peaks in the transmission spectrum 
centred at the frequencies around ±1^. Inhomogeneous 
broadening will damp the collective Rabi oscillation and 
lead to the line broadening of the transmission peaks. 
In Fig.3(a-b), we show that the linewidth of the trans- 
mission peaks is greatly narrowed with continuous driv- 
ing field, because the inhomogeneous broadening is sup- 
pressed. Once the effective collective coupling Ct is larger 
than the inhomogeneity of the driving field Bmax — Bmim 
the transmission function shows two peaks, the linewidth 
of which is mainly determined by the damping rate of the 
cavity and individual spins. Such a behaviour is very sim- 
ilar to the ideal rectangular spectral distribution [9]. To 
compare the effect of inhomogeneous broadening on the 
collective Rabi oscillation, we assume that the damping 
rate = and 7 = 0. In Fig. 3(c), it can be seen that the 
coherence time of the collective Rabi oscillation is greatly 
improved with continuous driving. 



V. LONG-LIVED QUANTUM MEMORY 
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Figure 4. (a) The state fidelity F at time t = 50 as a func- 
tion of the cavity-spin detuning S. (b) The time evolution of 
the state fidelity F(t) for a Lorentzian quantum memory with- 
out and with continuous driving. The detuning 6 is chosen as 
the optimal value. 



plains why a Lorentzian distribution of spin ensemble can 
not benefit much from the cavity protection, in particular 
as compared to the rectangular spectrum profile. 

The other critical factor for a long-lived quantum mem- 
ory is the damping of the cavity. To overcome this prob- 
lem, one can choose the cavity-spin tuning 6 such that 
the effect of cavity damping is suppressed but the mem- 
ory can still benefit from the cavity protection [8, 9]. The 
microwave photon can be stored to the dressed polariton 
mode as 



l^o(^)) 



(16) 



where cot (^) = {6/2Q) and \S) = {l/n)Ek9k\W' 
Thus, if the detuning is large the excitation is more likely 
to stay in the spin superradiant mode and suffers less 
from the damping of the cavity. While the cavity pro- 
tection against the inhomogeneous broadening requires a 
small detuning so that the effective collective coupling is 
strong. The storage time of a quantum memory achieves 
the best value for an optimal detuning [9]. The perfor- 
mance of a quantum memory is characterized by the state 
fidelity F{t) = |/(t)p with fit) = {M^)\(t){6,t)), where 
the state at time t is \(j){d,t)) = exp (—itH) \(j)Q{S)) = 

Qfo |1, + ak |0, t/c)- To calculate the state fidelity, 

one needs to solve the Schrodinger equation [8, 9] for the 
coefficients a/e, namely 



The strong coupling between the cavity and the en- 
semble induces an energy gap between the dressed po- 
lariton modes, i.e. superpositions of the cavity mode ac 
and the ensemble superradiant mode b = ^]^{gk/^)(^k^ 
which are protected from the dephasing caused by the 
inhomogeneous broadening. With such an effect of cav- 
ity protection, it is possible to build a quantum memory 
with a storage time much longer than the lifetime of the 
cavity mode and the dephasing time of the spin ensem- 
ble. The condition for the cavity protection is that the 
spectral distribution should decay sufficiently fast so that 
all the spin transition frequencies are far from the energy 
gap (arising from the collective coupling) [8, 9]. This ex- 



<^o(^) = -ioJcC^oit) - i^gk(y-k{t) (17) 

k=l 

c^kit) = -igk(^o{t) - iuJk(^k{t)' (18) 

where Qc = oOc — uok = cj/c — ^7/2. By performing 

the Laplace transformation, we obtain 

N 

sAo(5) - Ao(0) = -iuJcAq^s) - i"^gkAk{s) (19) 

/c=l 

sAk{s) - Afe(O) = -igkA^i^s) - zcDfeAfc(s). (20) 
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The solutions of ^0(5) and Ai^{s) are 



Aois) = 
Akis) = 



cos(-) + ^sin(-)M(s) 



T{s) 
9k 



■ ii^k 



— sin(-)- 

n ^2' s + iLOk 



(21) 
(22) 



where 



N 



S -\- lUJj 
= 1 



s^iCjc^ M{s) ' 
(23) 

The Laplace transformation of f{t) = {(j)o{S)\ (j){S^t)) = 
cos(f)ao(t) - {iMMDEkdkMt) is 



AT 



fis) = cos(-) Ao(s) - (l/n) sin(-) ^ gMs) (24) 



k=l 



from which one can calculate the state fidelity F{t) and 
compare the storage time of a quantum memory without 
and with continuous driving. 

We have compared the state fidelity for a Lorentzian 
distribution of spin ensemble with and without contin- 
uous driving, namely calculated the solution in Eq.(24) 
for the original Hamiltonian in Eq.(13) and the dressed 
Hamiltonian in Eq.(14). It can be seen from Fig.4(a-b) 
that: (1) The increasing collective coupling has very lit- 
tle effect on the state fidelity for Lorentzian distribution 
if there is no driving as expected; (2) The storage time 
can be improved by orders of magnitude with continu- 
ous driving; (3) The continuous driving also enables the 
quantum memory to benefit from the cavity protection, 
which suppresses the effect of the inhomogeneity in the 
driving field and thereby improves the performance with 
increasing collective coupling strength. We stress that 
the parameters in our numerical simulations are within 
the state-of-art cavity QED technologies. Careful opti- 
mization may easily improve homogeneity and strengthen 
the drive. To be more specific, we conservatively take 
the inhomogeneous broadening as A = IMHz and the 
damping rate of the cavity k = O.lMHz (i.e. the cav- 
ity lifetime is 10 jas). The storage time of a Lorentzian 
quantum memory is limited by the dephasing of the 
spin ensemble (about l/is). With continuous driving 



Bmin = 10(20)MHz and B^ax = 10.5(20.5)MHz, the 
state fidelity can be as high as 0.5(0.7) even after 50/i5. 



VI. SUMMARY AND DISCUSSION 

In summary, we have shown that continuous driving 
can suppress the effect of inhomogeneous broadening on a 
quantum memory based on spin ensembles. Besides, the 
obvious role of decoupling individual spins from magnetic 
noise, continuous driving can also reduce and reshape the 
spectral distribution. Even if the continuous driving field 
itself is inhomogeneous over the ensemble, it can make 
the tails of the spectral density profile decay much faster 
and approach the ideal rectangle distribution. This fea- 
ture then allows us to take advantage of the cavity pro- 
tection in the strong coupling regime and to achieve much 
improved performance of a quantum memory. We con- 
sider the Lorentzian spectral distribution as an example, 
as it is relevant in many solid-state systems, e.g. spin 
ensemble in diamond. Our results demonstrate that the 
linewidth of the double peaks in the transmission spec- 
trum is greatly narrowed by continuous driving. Further- 
more, the storage time of a driven quantum memory can 
be extended by orders of magnitude as compared with 
the one without driving. We stress that the main feature 
of our scheme is that the effect of the inhomogeneity of 
continuous driving field can be suppressed by the strong 
collective coupling via cavity protection. As comparison, 
one could invert all the spins half way through the stor- 
age period to eliminate the inhomogeneous broadening, 
nevertheless, the pulse inhomogeneity will inevitably lead 
to inversion errors over individual spins which can not di- 
rectly be compensated by cavity protection as occurs for 
continuous driving. We expect that the idea of driven 
quantum memory can be extended to other physical sys- 
tems, e.g. ion Coulomb crystals [18, 19]. 
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